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Abstract 

This paper derives the MKN (Maung-Kahana-Norbury) theory of bound states which 
incorporates the Lande subtraction technique to remove the singularities of the Cornell 
potential. 
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1 NRSE in Momentum Space 



Non-relativistic Schrodinger equation (NRSE) in configuration space has been solved ex- 
actly for some potentials, such as the Coulomb and simple harmonic oscilator potentials. 
NRSE with a linear potential can be solved analytically for the S'-state only as we will 
show later. For / > 0, we resort to numerical methods. NRSE r-space codes are com- 
monly known to be conditionally unstable [|I], 0], while the momentum space codes do not 
have the same problem. The momentum space code has an additional advantage of being 
easily adaptable to relativistic equations. NRSE in momentum space takes the form 



2jl 



y(q)0(p')rfp' = E0(p) 



where q = p — p'. 
Proof: 

NRSE in momentum space can be derived from its configuration space counterpart 



by Fourier transform. First we define the following: 
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and ignore factors of 2tt in Fourier and inverse Fourier transforms. As usual, we assume 
periodic boundary condition or V0 = = at infinity. We Fourier-transform Eq. 
term by term. The term on the right hand side of Eq. ^ is obtained simply by Eq. 
The first term involves Sl'^ip and is transformed as 
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The second term involves V"(x)-?/'(x). With Eq. Q], the second term of Eq. 
formed as 



is trans- 



/ 



K(x) 



//V(x)„p>.<---.3.,V 
j V{p - p')^{p') d^p' 
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At last, Eqs. 0, and are all needed to Fourier transform Eq. into Eq. |[T[]. The 
proof is complete. 

In this paper, we attempt to re-derive the results previously obtained by Maung et al. 
0] in 1993. The power law potential in r space is given by 
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Let G = h = c = 1. Define q = p — p'. The momentum space potential can be 
obtained by Fourier transform. 
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The final form of the momentum space potential is 
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where = — 1 corresponds to the Coulomb potential and = 1 the linear potential. 
Together they give the Cornell potential 



V{q) = V^^(q) + V^^(q) = V^=-^(q) + V^^='(q). (16) 

Next we want to perform a partial wave expansion of . There are 3 useful formulas, 
the Wigner-Eckart Theorem Q 

< E'l'm'\T\Elm >= 5n5n.'ME), (17) 

the addition of spherical harmonics 

9/ -I- 1 

Y.^U^)^LiSi') = ^^'(^o^^)' (18) 

m 

and the orthogonality of spherical harmonics 

dQY;^{Q)Y,^iQ) = 6in6m'„., (19) 
which are used in deriving the following result. 
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= XI 5Z < ^1 < ^1^^ lm\V^\l'm' >< l'm\Q' > \p' > (22) 

Im I'm' 

= X < f]|/m ><p|V^^|p' >< /m|fi' > (23) 
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Im 



9/ _|_ 1 

= E^^^''(^''^'')^^(cose) (25) 
I 

In scattering and bound state problems, it is customary to expand the momentum space 
wavefunction 0(p) in partial waves, such that 

0(P) = ^Cnlm4>nl{p)yim{^), (26) 
nlm 

where Cnim^s are coefficients of the expansion [p. 396 of Ref 1 ]. 



The non-relativistic Schrodinger equation in momentum space is given as 



^ ~ C j 0(P) = y l^^(q)0(p')rfV (27) 
= J <p\V^\p' >(p{p')d^p' (28) 

Expand NRSE in partial waves. 

^-f^)$^C„,^0nKp)>^^rn(f^) = j v'^ dp' d^' Y.V!' {v, V)YUmLm 
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The nZ-th terms can be separated by inspection with the help of the identity E (pniip) = 
Eni 4>ni{p)- The partial wave NRSE is 

(Eni - Mp) = J p"dp'Vi''{p,p')Mp)- (31) 
Use the orthogonality of Legendre polynomials, 

j Pi{x)Pv{x)dx = i^^^i'h (32) 

and Eq. (^), we can calculate the potential matrix elements as follow. 

j <p\V^\p' > Pi{cose)dcose (33) 

2/' + 1 /"^ 
= y — —Vi^{p,p') / Piicos 9) Pi> (cos 9) d cose (34) 

= ^^''(P,P') (36) 

In order words, 

Vi^{p,p') = 2n J V^{ci)Pi{cos9)dcos9. (37) 



Define 



y = 



2p'p 



and use the definition of the Legendre polynomial of the second kind Qn{z) 
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we can modify Eq. ( pT]) as 
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The coulomb potential corresponds to iV = — 1 and has the form 

TT V^O p'p 

The linear potential corresponds to = 1 and has the form 
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There are 3 useful relations in terms of Legendre polynomials of the 2nd kind 
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The singularities of Vf^lpjp') and Vi^{p,p') come from the singularies of Qi{y) and 
Q'l'iy). From Eq. (^0]), it is obvious that the singularities of Qi{y) and Q'liy) again come 
from those of Qo{y), Q'o{y) and Qoiy). In order to treat the singularities of the momentum 
space Cornell potential, we need to control the singularities of Qo{y), Q'oiy) and Q'q{ii) 
first and foremost. Substitute Eq. ( ^81) into Eq. (|^), we have 
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Differentiating Eq. 
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Differentiating again gives 
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There are two useful identities which we want to prove: 
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Proof: 

The integral in Eq. is derived as follow: 
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The results of Eq. ||6^ come from relations BI((256))(10) and BI((256))(11) in Gradsliteyn 
and Ryzhik 0. 

Eq. ||56| has 4 simple poles: a = p + ^r/, a*, /5 = —p + ir] and P*. Qo(y) can be rewritten 

as 

■1 1 
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The contour integral f QQ{y)dz over the upper complex plane has 2 residues: Res{a) and 
Res{P). Use the formula Res{zo) = lim2^^„(2; — zo)f{z) to calculate these residues, 
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Since the contour at infinity is zero and QQ{y) along the real axis is symmetric around 
the origin, we obtain 



Qoiy)dp' = ^(f^Q'o{y)dz 
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[rf' / p' p)Qq{ii) has the same poles as Qo(i/) but of order 2. Residues of order m are 
calculated by the formula 
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Again we can simplify the algebra by rewriting Eq. |^ 
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Hence the sum of residues is 
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Since {tj^ /p'p)Qo{y) is symmetric around the origin, we can integrate along the same 
contour as before and obtain 
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From these resuhs, it is obvious that Eq. |^ is true. The proof is complete. 

A simple example is the momentum space Schrodinger equation with a linear potential 
in the S'-state ^, 



2/i Trp^ Jo 



(f>no{p') dp' = EnO </)„o(p), 
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v^(p,p') 



where y = {jP' + p'"^) I2p'p. Lande subtraction ^ involves subtracting a zero term 
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from Eq. |^ such that 
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Using Eqs. , the integral in Eq. | |100| ] for p > in the limit of ?/ — 1 can be shown 
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The order of the hmits in Eq. [p.01|| is important. The reverse order will lead to the 
nonsensical result / Qoiv) dp' = 0. Next, in the limit of p,p' — >■ 0, (p + p'Y = (p — pY- 
By substituting this equality into Eqs. | |5^J55|] , it can be shown again that the integral in 
Eq. ||100|| vanishes for p ^ aX y = 1. At the end, the integral vanishes at y = 1, Vp. 
Away from the singularities, both integrands in the integral of Eq. [[100|| are finite. By 
taking rj —>■ 0, the first integrand vanishes. The final form of Eq. [|10CI|| is 



Q'oiy) [0no(p') - 0no(p)] dp' = E„o0no(p), 



(102) 



where Q'^iy) = l/{l-y^). 

The momentum space NRSE with a coulomb potential is given as 
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Use Eq. |Q to subtract out the logarithmic singularity and obtain 
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Before we perform Lande subtraction on the NRSE with a confining potential, we need 
the identity 

P/(l) = (105) 

Proof: 

We use the recursion relation xP^ — P[_^{x) = lPi{x), the equality Pi{l) = 1 to obtain the 
following relations: 

pui)-pL2a) = 1-1 



piii 



P[{1) = 2 
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Add these relations and use Po{x) = 1 or Pq{x) = 0, we prove Eq. [|105|| . Next we want 
to examine the singularities of Q'i{y) and Q'l{y). Differentiating Eq. ||5^ once and twice, 
we have 
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rfPiQa, rfPlQ'Q and 'rfwi_^ vanish in the hmit of — > 0. The momentum space confining 
potential is 
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Pl{y)Q,{y) - w[_,{y) 



(109) 



Lastly, perform Lande subtraction and use Eq. [|105|| . Take the limit of = 0, we derive 
the momentum space NRSE with a confining potential as 



„2 \ roo 

^Mp) + / Pi{y)Q',{y) Mp' 



iP) 



dp' 



+ 



A, 



p 



T^P" Jo 

^ — ^ — itpMp) 



P'<Pnl{.P') 



Piiv) 

l{l + 1) p(f)nl{p) 

2 Pliy) 



dp' 



X, 



-^P Jo 



w'i_MMp') dp' = E^iMp). (110) 



2 Exact Solution of NRSE with a Linear Potential 



In the next section, we are going to solve the NRSE numerically with a linear potential. 
In this section, we will solve the same equation exactly so that we can use the analytic 
results to check their numerical counterparts. The Hamiltonian equation can be written 

as 

' rf2 2 d\ 

+ _ R-2fi[XLr-E]R = 0. (Ill) 
^ar^ r dr J 

Let S = r R, then Eq. ||111|| can be simplified as 



d^ 

—S-2fi[XLr-E]S = 0. 
dr^ 



If we define as new variable 



I ^yiX^r-E], 



Eq. ||113|| can be transformed as 



S" — xS = 0, 



(112) 

(113) 
(114) 



which is the Airy equation. The solution which satisfies the boundary condition of 5 ^ 
as X — * oo is the Airy function Ai(x). Fig. |jl| illustrates the graph of the Airy function. 
By noticing that S = rR vanishes at r = 0, we infer that S{r = 0) must coincide with 
a zero of Ai(x). It is made possible by letting the eigen-energy act as a horizontal shift 
which shifts the origin to the left along the x-axis. If S is plotted against r instead of 
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X, S will vanish at the origin if En is chosen appropriately. This conclusion leads to the 
eigen-energy formula 



where x„ is the n-th zero of the Airy function counting from a; = along —x. 

In Norbury et. a/.'s ^ paper, the values Al = 5 and ^ = 0.75 are used. The eigen- 
energy formula is 



Table | hsts the zeros of the Airy function and the corresponding exact eigen-energies. 

3 Conclusion 

The p-space formalism shown in this paper can be applied to any arbitrary potential in 
principle although only the cases of n = — 1, 1 in r" are considered in this paper. More 
complicated potentials require the calculation of integrals involving other powers of r. 
The potential involving the coupling constant a{r) in QCD near the asymptotic freedom 
region is such an example. In these cases, we may need to reply on numerical integration 
to evaluate the integrals of the Lande subtraction terms. 
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Tabic 1: Zeros of the Airy function and the corresponding eigen-energies in GeV for I — 0, 
PL = 0.75 GeV, Al = 5 GeV. 



n 




En 


1 


-2.33810741 


5.972379202 


2 


-4.08794944 


10.44211404 


3 


-5.52055983 


14.10152355 


4 


-6.78670809 


17.33572806 


5 


-7.94413359 


20.29221499 


6 


-9.02265085 


23.04714148 


7 


-10.04017434 


25.64626764 


8 


-11.00852430 


28.11978666 


9 


-11.93601556 


30.48893766 


10 


-12.82877675 


32.76937540 
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Figure 1: The graph of the Airy function-The zeros are all negative. Since S = r R, S 
must be zero at the origin r = 0. The eigen-energy acts like a horizontal shift which shifts 
the origin back to a zero. 
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